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Abstract 

The cosmic microwave anisotropics in a scenario of large scale structure formation 
with cold dark matter and texture are discussed and compared with recent observa- 
tional results of the COBE satellite. A couple of important statistical parameters are 
determined. The fluctuations are slightly non gaussian. The quadrupole anisotropy is 
1.5±1.2 x 10 -5 and the fluctuations on a angular scale of 10 degrees are (3.8±2.6) x 10~ 5 . 
The COBE are within about one standard deviation of the typical texture + CDM model 
discussed in this paper. Furthermore, we calculate fluctuations on intermediate scales 
(about 2 degrees) with the result AT/T(6 ~ 2°) = 3.9 ± 0.8) x 10~ 5 . 

Collapsing textures are modeled by spherically symmetric field configurations. This 
leads to uncertainties of about a factor of 2. 



1 Introduction 



Global texture is the latest of a set of models based on a simple physical idea: the universe 
begins in a hot, homogeneous state and then, as it cools, undergoes a symmetry breaking 
phase transition that leads to the formation of topological defects [§]. These defects induce 
perturbations that seed the formation of galaxy and large-scale structure. Unlike stable 
defects (cosmic strings, domain walls and monopoles), textures are short-lived. In a light- 
crossing time, they collapse producing a singularity for one instant of time and radiate their 
energy in a spray of Goldstone bosons. They do, however, persist long enough to leave an 
imprint as compensated perturbations in the matter and radiation density. The simplicity 
of the texture scenario is appealing, the model has only a single adjustable parameter, 
the scale of symmetry breaking, which we normalize to obtain the observed galaxy-galaxy 
correlation function. 

In a series of papers, we and our collaborators have examined the texture-seeded SI = 1 
CDM-dominated cosmogony and found promising results. Analytical studies of galaxy 
formation in this scenario suggest that the model has several highly attractive features ||: 
early star formation, which could reionize the universe and smear out fluctuations in the 
microwave background; early formation of galactic spheroids could account for their high 
mean density without significant dissipation and would provide a fertile environment for the 
formation of a massive black hole that could power quasars; a galaxy mass function whose 
slope and amplitude is consistent with observations; and galaxy formation at moderate 
z. Numerical simulations of the growth of large scale structure in the texture scenario 
find striking consistency between theoretical predictions and observations of large scale 
structure Q, [q, g: the galaxy-galaxy correlation function is consistent with observation; 
the derived galaxy mass function is well-fit by observations; the model produces coherent 
velocity fields on large scales and the angular correlation function of galaxies is consistent 
with observations. 

In this paper, we calculate the signature of microwave background fluctuations induced 
by texture on the scales probed by COBE. In section 2, we present our formalism for 
calculating the response of dark matter, baryons and photons to the collapse of a single 
spherically symmetric texture in an expanding universe. These calculations extend earlier 
work ^] that calculated the fluctuations induced by a texture in flat space. In section 3, 
we describe how we sum over the contributions from many textures and make a synthetic 
map of the microwave background. Section 4 presents our numerical results and conclusions. 

2 The Effects of the Collapse of a Texture 

Texture arises when the symmetry breaking "Higgs" field chooses to lie at different places 
on its degenerate vacuum manifold in different regions of space during the phase transition. 
This can produce configurations in which the Higgs field has non-trivial winding number 
on the vacuum manifold, which we call "knot". Such a knot is unstable to collapse — it 
shrinks until the gradient of the Higgs field becomes large enough that the field lifts off the 
vacuum manifold for an instant, undoing the knot. This happens when the knot shrinks to 
a microscopic size, of order the inverse symmetry breaking scale. After the knot unwinds, 
the field is free to radiate away completely. 

The cosmological evolution of texture is remarkably simple: in any given horizon volume, 
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there is a fixed probability of the field configuration being wound into a knot [[I], M. Once 
there is sufficient winding number density within the horizon, the knots collapse at the 
speed of light. As the horizon size grows and "connects" new regions, the process repeats in 
a self similar manner, with new knots continually collapsing on progressively larger scales. 
As they collapse, unwind, and radiate away, the knots generate a gravitational field that 
perturbs the surrounding matter. 

In this paper we study the simplest case of texture: that formed by a completely broken 
global SU(2) symmetry. It is important to notice, however, that the general picture of 
cosmological texture evolution just discussed is expected to be the same regardless of the 
"species" of texture. Thus we expect many of the main features of the resulting structure 
formation to be species independent. Furthermore, it turns out that for any texture the 
dynamics of the field is dependent only on the geometry of the vacuum manifold, and cannot 
be affected by the shape of the Higgs potential or even the symmetry breaking scale. The 
latter does enter our scenario as the one "tunable" parameter, but affects only the amplitude 
of the induced fluctuations, not the pattern. We will see that normalizing the fluctuations 
such that the resulting cosmic structure matches well with observation implies a symmetry 
breaking scale quite natural for GUT textures. 

2.1 Formalism 

In this section, we will present a formalism for calculating the response of matter and 
radiation to the collapse of a single texture knot. Throughout this paper, we will assume a 
flat Friedman universe, Qtot = 1> dominated by CDM and we work with conformal time so 
that the background metric is given by 

ds 2 = a 2 (-dt 2 + dx 2 ) . 

We will model the effects of texture by considering the response of matter to the collapse 
of a single spherical texture. This approach has proven fruitful in describing the response 
of dark matter to textureQ. The spherically symmetric ansatz for an S 3 texture (i.e. a 
configuration of a scalar field living on S 3 which winds once around S 3 as r goes from to 
oo) unwinding at a given time t = t c is ]?]] 

<fi = r/ 2 (sin xsin# cos <p, sin xsin# sin if, sin x cos 9, cos x) , (1) 

where 6 , (p are the usual polar angles and x( r > t) has the following properties 

X {r = 0,t<t c ) =0 
X(r = 0,t>t c ) =vr 
x(r = oo,t) =7T. 

Since the texture seed is already a perturbation, it can be evolved according to the back- 
ground equation of motion, □</> — ((f)-0(f))cj) = 0, which yields 

d 2 tX + 2{a/a)d tX - d 2 rX - -d rX = • (2) 

The energy momentum tensor of the texture field, T^ v = d^4> • dv<t> ~ \9tivd\4> " d x <p, can 
be parametrized in terms of dimensionless functions, 

Too = Wfp/l 2 , Toi = -Wfv,i/l , 
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where 



f p /l 2 = l/8[(d tX ) 2 + {drX? + ^r^] , (3) 



f P /i> = i/mxf - \{drx? - 2j ^\ , (4) 



Ull = \\ (d tX )(drX)dr , (o) 
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A , ! r / P n2 si^Xi . 3 r 2 sin 2 *^ 

^U = -r[(drX) — +7/ l( d rX) — ] — ' (6) 

4 r z 4 Joo r r 

The length parameter Z is introduced to keep the functions f p to f n dimensionless. A 
convenient choice is to set I equal the typical size of the perturbations, e.g., the horizonsize 
at texture collapse. f p /l 2 and f p /l 2 are the energy density and isotropic pressure of the 
texture field. f v is the potential of the velocity field and f n is a potential for anisotropic 
stresses. 

Cosmological perturbations are usually split into scalar, vector and tensor contributions 
which do not couple to first order (here the terms scalar vector and tensor refer to the 
transformation properties on hypersurfaces of constant time). The split into scalar, vector 
and tensor perturbations is nonlocal and thus acausal. But spherically symmetric pertur- 
bations are of course always of scalar type (due to the adoption of spherical symmetry, we 
loose, e.g., all information about gravity waves produced during the collapse.). In this case 
it is very convenient to apply gauge-invariant cosmological perturbation theory which is 



adapted to a split into scalar, vector and tensor modes [1C, 11, 12, H|. In this paper we use 
the formalism and notation introduced in ||. 

The gravitational field can be described in terms of the gauge-invariant Bardeen po- 
tentials and ^. The explicit definition of the Bardeen potentials in terms of the metric 



perturbations is given in [10, 11, 12, 15]. For perturbations which are much smaller than 
the size of the horizon, I <C Ih, ^ corresponds to the Newtonian potential, and $ is related 
to the perturbation of the 3-curvature on hypersurfaces of constant time: 

-(4/a 2 )A$ = 6R® + 0(l/l H ) ■ 

Hence for sub horizon perturbations of ordinary matter (|T 00 | 3> |T l:, |), the well known 
Newtonian limit of Einstein's equations (see, e.g, |l3|]) yields 

A$ = -A* . 

Einstein's equations determine the Bardeen potentials induced by the matter perturba- 
tions H . In addition to texture we want to describe dark matter and radiation. Since dark 
matter has zero pressure, we just need a variable describing its density perturbation. 

By Pd we denote the density of dark matter. D is a gauge invariant combination of the 
density perturbation, 5 = (pd — Pdo)/ Pdo, the potential for the peculiar velocity, Iv , Vi = 
—lv,i and geometric terms. 
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In this variable, the evolution of dark matter fluctuations is governed by 

D + (a/a)D - ^G Pd a 2 D = e(f p + 3f p )/l 2 , (7) 

where e = l&irGrj 2 . A derivation of this equation is presented in||. After decoupling, 
the evolution of radiation is given by the one particle Liouville equation. Its first order 
perturbation yields a differential equation for A4, which is a gauge-invariant version of the 
perturbation of the energy integrated photon distribution function. 



47T 



M = — Tp 6 dp . (8) 
p Jo 



M evolves in response to fluctuations in both the $ and Vl/ potentials: 14, 15] 

M + J%M = 4rf%(^ - *) . (9) 

The unit vector -y denotes the direction cosines of the photon momentum. In spherically 
symmetric systems, the photon evolution equation @ reduces to, 

1 - u 2 

M + pd r M + — d u M = Apdri® - m (10) 

r 

where [i is the direction cosine of the photon momentum in the direction of r. 

The photon evolution equation is more transparent in characteristic coordinates, (t, r, b) 
where b = ryl — p? is the impact parameter and r = t — rp is the time at which the 
photon is at its minimum distance from the texture, the 'impact time'. In these variables 
fllC|) simplifies to 



d t M(t,T,b) = Apd r {>& - , (11) 

so, that 

M(t,r,b) = M{t ,r,b) +Ap C dt'd r (®-y) . (12) 

To write down the perturbed Einstein equations we have in principle to calculate the 
energy momentum tensor of radiation from M.. But since we are only interested in late 
times, pd 3> j0 7 , we may neglect the contribution of radiation to the density perturbation. 
The potential is then determined by the texture and dark matter perturbations alone ||: 

A$ = -e(/ p // 2 + 3(a/a)f v /l) - ^Ga 2 p d D , (13) 

On the other hand, since dark matter does not give rise to anisotropic stresses, we have 
to take into account the contribution of radiation to the latter. To calculate the anisotropic 
stresses of radiation we recall the definition of its amplitude, II (see e.g. [14, [0]]): 

STi - ^ST}S{ = P [n;l - i AUSf] . 

The anisotropic contributions to the energy momentum perturbations of photons are given 
by 
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Using these equations and spherical symmetry one finds 

II" - n'/r = | fy - l -)Md^ = : ^M 2 (r, t) . (14) 

This anisotropy and the anisotropy of the texture field contribute to the sum of the two 
Bardeen potentials || 

A(^ + $) = -2A(eU + ^7TGa 2 Pj U) . (15) 

Inserting the results for d r (& — x I / )(r(r, fx), t) at a given time step, we can evolve the 
brightness function A4 according to (|l2|). The function Ai is very simply related to the 
temperature perturbation in a given direction in the sky: To see this, let us express the 
perturbation of the photon distribution function as a perturbation of the temperature: 
/ = /(o) + $f = f(0)(jP/T(t, r, 7)) so that 5f = ^^-5T/T. If we know neglect the difference 
between gauge-invariant and gauge dependent variables and set 5f = T ', we find from (||) 
the simple relation Ai = 45T/T. In the gauge-invariant framework, this just means that 
.M/4 is a gauge-invariant variable describing the perturbation of the temperature in a given 
direction (for a more explicit discussion see |15|| ). 

We would like to choose initial conditions that are physically plausible. If the phase 
transition that produced texture occurred in an initially uniform universe, causality requires 



that there are no geometry fluctuations well outside the horizon [16]. This implies that 
initially, at t, <C t co i, we must require ^ = <I> = 0. We want to compensate as much as 
possible of the initial texture fluctuations with an initial dark matter perturbation. Hence, 
we use as initial conditions for the density field, 

D(r,t = U) = - 47T G a 2 pd (fp/ 12 + H^Mfv/l) ■ (16) 

This initial condition implies that metric fluctuations are induced by the differences between 
the texture equation of state and the equation of state of the background matter. This choice 
yields $ = at t = ti, but not ^ = 0. Due to its equation of state the dark matter cannot 
compensate the anisotropic stresses of the texture. We thus must compensate them by an 
initial photon perturbation. For \I/ to vanish, we have to require according to (|l~5|) 

n = *— fw . 

4-7rGa 2 p 7 

Clearly this does not lead to a unique initial condition for Ai, but if we in addition require 
the zeroth and first moments of Ai to vanish it is reasonable to set 

M ^*-^—5^y-\^-\^- <i7) 

Together with initial conditions for the texture field x( r >^ = t ne requirements (|l6| ) 
and ( |i~7| ) and the evolution equations (|2|) to (^) determine the system which we solved 
numerically. 
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2.2 Texture Collapse in Flat Space 

Before we present numerical methods and results, let us briefly discuss the flat space case, 
where, if we ignore the response of dark matter to the texture, the microwave fluctuations 
can be calculated analytically. The Bardeen potentials of a self similar texture collapsing 
at t c = are |§ 

e, ,r 2 + t 2 , e, ,r 2 + t 2 x 



4 v I 2 ' ' 4 v t 2 ' ' 

so that (0) turns into 

d t M = w ^^ T2 =g{t^b). (18) 

Integrating that over time one finds 

r It — t It- — t 

M(t,T,b) = Ai(tj, r, b) + 2e[ , =(arctan , = — arctan ■ * ) 

K J { ' Vr 2 + 26 2V Vr 2 + 26 2 Vt 2 + 26 2; 

-lo g ( (2 ^ T)2+T2 + 262 )1 (19) 
gl (2t i -r) 2 + r 2 + 26 2jJ 1 ' 

a We are interested in the change of the distribution function due to the collapsing texture 
and thus want to perform the limit t — > oo , U — > —oo, with lim^—j—oo Ai(ti, r, b) = 0. The 
log term in ([II]) then contributes a constant C = lim(log(t 2 /t 2 )) which crucially depends 
upon how we perform this limit. It can take any value — oo < C < oo. This can happen 
because the flat space texture is an unphysical infinite energy solution. If we would change 
the energy momentum tensor in a way that the texture would be 'born' some time in the 
finite past, or if we would compensate it in a consistent way as we do it in the expanding 
universe, this problem would disappear. For simplicity, and since we are not interested in a 
r and b independent constant, we just perform the limit t = —ti — > oo which yields C = 0. 
Then the final change of the distribution function is 

^• t) = 2 ™77kw (20) 

Since M describes the fluctuations in the energy density, AT/T = AAi/4. 

In a flat universe the signal from a texture collapsing at t = as seen from an observer 
at time t and distance r Q would thus be 

= «r/g . U - r '^ 6 ■ (21) 
1 yj (t - r Q cos 6) 2 + 2r Q sin 2 6 

where 9 is the angle between the direction of the texture and the line of sight. A similar 
calculation using a specific gauge but performing a somewhat unphysical limit is presented 



in [17]. Unlike in the expanding universe (see 2.3), there is no horizon present in these 
calculations, thus, photons that pass the texture long before or after collapse, \t Q \ ^> r Q are 
still influenced by it and yield even a maximum temperature shift, 

AT 

— = ±evr/2 . 

In the expanding universe we expect AT/T to achieve a maximum for r Q ~ t Q — t c and to 
vanish for t Q — t c r Q . 
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2.3 Texture Collapse in the Expanding Universe 

The coupled set of equations that describe the evolution of texture, the fluctuations in 
the dark matter, the fluctuations in the metric, and the fluctuations in the microwave 
background in the expanding universe must be integrated numerically. We use a second- 
order accurate leapfrog scheme to advance these equations. 

We will assume a specific set of initial conditions for the texture field. One of the 
sources of uncertainty in our calculations of the microwave background is the choice of initial 
conditions. We have found that different choices of initial conditions can yield microwave 
fluctuations that differ by ~ 50% from the values reported here. We choose initial values for 
X that confine the fluctuations in x to a region comparable to the horizon size at collapse: 

x = Ar/r + B(r/r ) 2 + C(r/r ) 3 r < ar (22) 

X = vr r > aro 

where A, B and C are determined by the requirement that x~ ^ an d its first and second 
derivatives vanish at r = arQ. For a = 4/3 this texture collapses at t ~ rQ. In contrary to 
the self similar flat space solution given in [f7|, ||], for this set of initial conditions, f n does not 
vanish initially and the ^(U) contribution must be compensated by a photon anisotropy II 
as explained above. 

Figure 1 shows the microwave background fluctuations induced by the collapse of the 
texture as a function of r for small impact parameter for both the flat space solution (dotted 
line) and the texture in the expanding universe. The dotted line in the figure is the flat space 
result from section 2.2. Figure 2 shows the microwave background fluctuations induced by 
texture collapse as a function of impact parameter. Note that the temperature fluctuations 
are induced only for photons that pass within the event horizon of the texture. 

2.4 Compensation 

In addition to the choice of sensible initial condition, we must keep in mind that spherical 
symmetry in an acausal concept and induces accelerations on super horizon scales. This is 
pictured by the dashed line of Fig. 3. There AT/T is shown as a function of the impact time 
r at constant time t ~ t c and impact parameter b ~ i c /50.. One sees that photons with 
impact time r ~ —40. = — t c , i.e. photons which had a distance r > t c from the texture at 
the time of the big bang, t = and always moved in the direction opposite to the texture 
(/j, ~ 1), still experience nearly half of the maximum energy shift even though they have 
never been in causal contact with the region of texture collapse (which at the collapse time, 
t = t c encompasses a radius of t c ). 

These acausalities would lead to an overestimate of AT/T. To prevent this we have 
to compensate in such a way that there are no accelerations on scales larger than the 
horizon. At each time we thus set the compensation radius equal to the size of the horizon, 
fcp = Ih = t. We then set f p (r) = f p (r) = f v (r) = /^(r) = for r > r cp . On small scales 
we subtract the mean of f p and f p , so that 

/ f P = f P = o ■ 

Jo Jo 

A texture compensated in this way does not produce any gravitational fields on scales larger 
than r co . 
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We have added this compensation to the calculation of the texture sources. The signal 
of a collapsing texture has then reduced to the solid line in Fig. 3. Here photons which have 
never been and will never be in causal contact with the texture (r < 0) and photons which 
are not yet in causal contact with the texture (r > 2t c = 80) do not acquire any appreciable 
energy shift. A result which is much more physical than the dashed line. 

Using this compensation scheme, we assume that spherical symmetry is a good approx- 
imation inside the horizon and that the stocastic texture field outside the horizon can be 
neglected. The question if this assumption is justified can only be answered by a full 3d 
simulation like |23|]. We believe that the issue of compensation and the arbitrariness in the 
initial conditions yield a factor of about 2 uncertainty in our results (which is much bigger 
than our numerical errors). It might well be, that by neglecting super horizon fluctuations 
completely, we slightly underestimate AT/T but most probably by less than a factor 2. 



3 Textures and the Microwave Sky 

In the previous section, we described how the collapse of a single texture produces fluc- 
tuations in the photon temperature. In this section, we sum the contribution of multiple 
textures and describe how to construct microwave maps of the night sky. 

Since COBE observations cover the entire celestial sphere, we construct a numerical grid 
consisting of 1 square degree patches. These patches are arranged so that they cover equal 
areas and each patch has roughly the same shape. 

Since the texture fluctuations are in the linear regime, we assume that the contributions 
of each texture to the fluctuations at each point in the sky can be added independently. 
We randomly throw down textures everywhere within the event horizon using the texture 
density distribution function Q, 

— (23) 
dt~2ht^ [Z6) 

follow the collapse of each texture and sum their contributions. We include only contri- 
bution of textures that collapse after recombination. Textures that collapse earlier do not 
contribute significantly to microwave fluctuations on scales accessible to COBE. 

In order to simulate the COBE observations, the map of the night sky is smoothed 
with a Gaussian beam with a FWHM of 7°, the angular resolution of the DMR detector on 



COBE[|T8, 18]. After computing and removing the intrinsic dipole contribution (10 -5 — 10 -4 ) 
and any monopole fluctuations, we then compute the microwave quadrupole, the r.m.s. 
pixel-pixel fluctuations and other statistics of the microwave sky. 
The simulated 'COBE-map' is shown in Fig 4. 



4 Results 

The amplitude of the microwave background fluctuations depends upon the scale of symme- 
try breaking associated with the texture. If the texture model is normalized so that it can 
reproduce the galaxy- galaxy correlation function, then e = 5.6 x lO -4 ^" 1 , where 6, the bias 
factor, is the ratio of the mass- mass correlation function to the galaxy- galaxy correlation 
function H, |24|. Comparison of the predictions of the texture model with observations of 
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clusters p5fl suggest that b « 2 — 4. Hydro dynamical simulations of texture-seeded galaxy 
formation|]5| suggest that b as 2. 

Figure 5 shows the distribution of quadrupole values calculated from 100 realizations of 
model. Averaging over the realizations, we find an r.m.s. value for Q of 

Q = (1.5 ± 1.2) x 10~ 5 e , 

where we have set eo = e/2.8 x 10~ 4 so that a bias factor of 2 and the above normalization 
would yield eo = 1. In 13 percent of the realizations the calculated value of the quadrupole is 
below the COBE result (Q = 5 x 10~ 6 ). Since only a handful of texture knots are the source 
of most of the large scale fluctuations, the quadrupole varies significantly from realization 
to realization. The 95% confidence range is 3.5 x 10- 6 e to 4 x KT 5 e . The COBE value 
is well within one standard deviation. 

The texture-induced microwave sky shows non-Gaussian features that distinguish it 
from the standard inflationary models. Figure 6 shows the distribution of pixel-to-pixel 
fluctuations calculated in one of the simulations. The dotted line shows a Gaussian whose 
dispersion matches the temperature dispersion of this simulation. The average temperature 
fluctuation averaged over 100 simulations is 

(AT/T) rms (6 = 10°) = (3.8 ± 2.6) x 10~ 5 e . 

This number is consistent with estimate based on three-dimensional simulations of texture 



induced microwave fluctuations [24, 23]. Figure 6 shows that the distribution of temperature 
fluctuations are only mildly non- Gaussian, the skewness of the distributions averaged over 
100 realizations is —4 ± 2.3 and the kurtosis of the distributions is 32 ± 29. 

Considering these results, the recently reported results from the COBE differential mi- 
crowave background radiometers [DMR] [JH]] are consistent with the texture + CDM sce- 
nario. The experiment has measured a value of Q = (0.5 ± 0.1) x 10 -5 for the microwave 
quadrupole and (AT/T) rms = (1.1 ± 0.18) x 10~ 5 for the pixel to pixel fluctuations in the 
microwave maps. Both these values lie within one standard deviations from the results of 
more than 100 realizations of the texture scenario. (Recent reexaminations of the COBE 
results suggest that the quadrupole may be somewhat lower f22[].) For the texture scenario, 
this implies that e ~ 3 x 10~ 4 and a bias factor of b ~ 2 are compatible with both, the 
large scale structure formation and the microwave background anisotropies produced in this 
scenario. 

Our simulation does not take into account the baryon/photon coupling before recombi- 
nation and does not properly model the process of recombination itself. It is thus not suited 
to account for textures which collapse before recombination, i.e., for angular scales beyond 
about drain = 2°. At this minimal resolution, the fluctuations in our simulations are 

(AT/T) rms (0 = 2°) = (3.9 ± 0.8) x 10" 5 eo . 

It is interesting that the cosmic variance (i.e. difference in values over our 100 realizations 
for this angle now has really dropped. This is clear, since more than 1000 textures have 
contributed to this signal in each map. Our number does contradict the south pole result 
|2C] but is in marginal agreement with the recent MAX experimant (third flight) plfl . In 



Fig 7 the autocorrelation function is shown with la upper and lower limits. On small scales, 
it is typically a factor 2-3 higher than the COBE result [IS]. This shows that the texture 



scenario (without reionization) marginally conflicts observations on small scales (9 < 5°). 
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In Fig. 8 the amplitudes of the spherical harmonics are presented. To specity our 
conventions, we define them here: 



(^/w^) = EE°'AW)' where 



(24) 



1=0 m=-l 



a lm = J dn(AT/T)Y^ . 



The amplitudes C\ are then given by 



Ci 



1 




(25) 



2Z + 2 



m=—l 



In Fig 8 we present in addition to the values < C\ > averaged over 100 realizations also 
their la upper and lower limits. Furthermore, the analytical fits for power law spectra with 
n = 1.5 (short dashed) and n = 1 (long dashed) are plotted. A simple investigation leads 
to a la limit of the spectral index of about n = 1.2 ± 0.3, although our numerical result 
represents definitely not a pure power law. 

Up to now we have only considered the statistical variations of our results which vary 
accidentally from realization to realization, i.e., cosmic variance. It is however important 
to note that due to the uncertainties in the initial configurations of the textures and the 
approximations inherent in modelling the texture as spherically symmetric, these estimates 
of the microwave background fluctuations are uncertain by a factor ~ 2, representing a 
systematic error. Since in this approximation we completely neglect all non topological 
fluctuations of the scalar field, we probably underestimate the induced fluctuations. On 
the other hand, a different choice for the scalar field leading to a different vacuum manifold 
M may as well alter the results by a factor of 2 (by altering the probability of texture 
formation). In addition, we throw down the textures randomly, uncorrelated whereas 3d 
numerical simulations hint that they are probably anticorrelated. 

Because of the non-Gaussian nature of the initial conditions, galaxy formation begins 
much earlier in the texture model than in Gaussian theories of structure formation. By 
z = 50, a fraction 6 x 10~ 4 of the mass of the universe has formed non-linear objects of 
mass greater than 1O 6 M — these objects may have reionized the universe. If 10% of the 
baryons in these objects are "burnt" in high mass stars, then the resultant energy, 60 eV 
per baryon is sufficient to reionize the universe. The optical depth of an ionized universe is 



so that we get an optical depth of 0.33 back to z = 50 in an fij, = 0.05, h 50 = 1, and 
£1 = 1 universe. A very rough estimate of the effects of reionization by smoothing each 
texture with a smoothing scale of about the horizon size at ^ ec , where z& ec is determined 
by T(z dec ) = 1, thus 



t(z) = 0.9- 10~ 3 ( 



0.05 



)[(l+z)3/'-l] 



(26) 
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This corresponds to an angular scale of 



Osmooth = t(z dec ) /t = (1 + z) 1 k, 5.7° 



( 



^50 \ 1/3 

0.05 ) 



(27) 



This smoothing length is somewhat larger than the resolution scale of COBE. If the forma- 
tion of objects leads to reionization prior to Zd ec this would suppress microwave background 
fluctuations on the smallest scales probed by COBE (~ 3°), but would not effect the fluc- 
tuations on the larger scales (> 10°) discussed above. 

We conclude again that our simulations indicate that the texture scenario can reproduce 
the CBR anisotropics measured by COBE and no pronounced inconsistencies with the 
COBE data are found. 

Due to the fact that a collapsing texture is a rare event, the standard deviations of 
fluctuations increase strongly towards larger scales. This leads to an increase by nearly a 
factor of 2 from the 1 — a lower limit of the fluctuation amplitude at 10 degrees to the 
corresponding limit at 2 degrees. This can be compensated by the damping of fluctuations 
due to reionization which can also affect scales which are somewhat larger than 6 smoo th- 

To further test this scenario it is important, in addition to a full three dimensional 
investigation of the CBR anisotropies, to discuss the damping by reionization in detail for 
anisotropies on intermediate scales, 1° to 5°. 

Our work confirms that the texture scenario is a viable candidate for large scale struc- 
ture formation which cannot be excluded by present observations. Its CBR anisotropies 
can be distinguished from those of a scenario with initial fluctuations induced due to an 
inflationary phase by two statistical criteria: The fluctuations are slightly non Gaussian and 
they increase somewhat towards smaller angular scales. A speciality of the texture scenario 
is furthermore the very distinct signal which one singular texture collapse would imprint 
on the microwave sky. Depending on the time of texture collapse and on our distance to it, 
we would see a roughly spherical ring of negative, positive or negative-positive fluctuation, 
see Fig. 9. Given the statistical distribution of textures (23), we expect to observe about 
5 — 10 such rings with an angular diameter around 10°. 

It is clear that although as a function of the impact parameter the phton signal shows 
a blueshift followed by a redshift, what is observed at a given moment in the sky is ei- 
ther dominated by redshifted or by bueshifted photons, since most photons arriving at the 
observer today have impact times which differ by much less than the typical size of the 
texture. Due to spherical symmetry, the signals form rings in the sky. In a more realistic 
3d simulation these rings would be deformed. The blueshifted rings are more affected by 
compensation and thus less pronounced as seen the map (Fig. 4). 
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APPENDIX 



A Damping by photon diffusion 



In this appendix we want to estimate the damping of CBR fluctuations in an ionized plasma. 
For simplicity we work here in physical, not confer mal time and we always consider a flat, 
matter dominated universe, so that t = to/(l + z) 3//2 and p-y < pb- 

For a perturbation with wave number k which is significantly smaller than the mean 
free path, the damping rate is then |E§] 



7 « k 2 t T /6 



(28) 



where tx = l/(o"Tn e ) is the mean free path of the photons. We now consider a texture which 
collapses at time t c = io/(l + z c ) 3 ^ 2 - The total damping, exp(— /), which this perturbation 
experiences is given by the integral 



/ 



-y(t)dt . 



(29) 



The end time t en d is the time, when the mean free path, tx equals t c , the size of the 
perturbation. We define Zd ec as the redshift when the photons and baryons decouple due to 
free streaming. This is about the time when the mean free path has grown up to the size 
of the horizon: txizdec) ~ i(^<fec)- To obtain exponential damping (t c < t en d), we thus need 
z c > Zdec- In this case damping is effective until 



1 Z e nd ~ 

and we obtain 



(1 + z dec )(l + z c 



f 



1 + Zdec 



l + Z c J 



3/2 



1 + Z C 
1 + Zde 



9/4 



(30) 



where we have set k 2 /6 = (2ir/t c ) 2 . In terms of angles this yields 

f(9) « 2{e/6 d f[(e d /efi 2 - 1] , 

where 6 d = + z dec ~ 5.7°. 

If the plasma ionizes at a redshift Zi, Zdec < Zi < z c , after the texture has already 
collapsed, damping is only effective after z\ if z en d < %i and instead of formula (^0|) we then 
obtain 



/ 



1 + 

l + z c 



3/2 



l + Z c 
1 + Z de 



9/4 



l±Ze\ 9/A 
1 + Zi 



(31) 



This damping factor can again be converted in an angular damping factor in the above, 
obvious way. The CBR signal of textures is thus exponentially damped only if 

Zi > Zdec ~ 107(0.05//i 50 ^b) 2/3 
and in this case only textures with 1 + Zd ec < 1 + z c < (1 + zi) 2 / (1 + Zd ec ) are affected. For 



12 



exponential damping after reionization at Z{ = 50 extreme values of the Hubble Constant 
and baryon densities far above that predicted from big bang nucleosynthesis and light 
element abundances are required. 

In this approximation we have neglected the the amount of damping which still may 
occur after z enc i- 
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Figure Captions 



Fig. 1: The hot spot — cold spot signal of a spherically symmetric collapsing texture in 
units of e ~ 2.8 x 10~ 4 . The horizontal variable r = t — r cos 9 denotes the 'impact time' (the 
time of closest encounter with the center of the texture) of a photon arriving at a distance 
r from the texture at time t traveling with an angle 9 with respect to the radial direction. 
The hot spot-cold spot is shown for photons with fixed impact parameter b = r sin 9 ~ 0.1i c 
(t c is the time of texture collapse) . 

A texture in a recombined expanding universe at t = t c , line (1) and t = 1M C , line (2) is 
compared with the flat space result (dashed curve). The second peak appearing at t = IM C 
is due to the dark matter potential. 

Fig. 2 The CMB perturbation in units of e ~ 2.8 x 10 -4 as a function of the impact 
parameter b for fixed r ~ 0.5i c = 10. The signal disappears at an impact parameter 
b ~ (1 — 1.5)i c (t c = 20 in the units chosen). 

Fig. 3 A comparison of AT/T from a compensated (solid line) and a not compensated 
(dashed line) texture is shown. The not compensated signal is acausal: a photon which 
had a distance « 0.5i c from the texture at the time of the big bang (t = 0) and aways 
was moving away from it (r « t c /2), i.e. never was in causal contact with the texture still 
experiences about half the maximum energy shift. These acausalities are removed in the 
compensated case. 

Fig. 4 A simulated COBE map as it might look in a scenario with texture + CDM. The 
color scheme goes from —4 x 10~ 4 (dark blue) to 4 x 10 -4 (deep red). The zero point is 
not the middle of the color palette but at the transition green/yellow. Monopole and dipole 
contributions are subtracted in this map. 

Fig. 5 The distribution of quadrupole values of 100 realizations of the microwave sky 
is shown in linear (a) and logarithmic (b) scale. The distribution is very broad with an 
average of 1.5 x 10 -5 and a standard deviation of 1.2 x 10~ 5 . 

Fig. 6 The statistical distribution of the microwave anisotropies for one realization of the 
microwave sky. The number of pixels with a given AT/T are counted. The distribution is 
slightly non Gaussian, negatively skewed (skewness = —3, curtosis = —1.2). 

Fig. 7 The auto correlation function is presented in units (/J.K) 2 (for a bias factor of 2, 
eo = 1, see text). In addition to the average values over 100 realizations (solid line), we 
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indicate the la upper and lower limits (dashed lines). 

Fig. 8 The amplitudes of the spherical harmonics for / = 2 to 64, averaged over 100 
simulations are shown (fat solid line) and the la upper and lower limits are indicated (thin 
solid lines). For comparicon we have also plotted the analytical fits for power law spectra 
with n = 1.5 (short dashed line) and n = 1 (long dashed line). 

Fig. 9 The spherical ring of positive-negative fluctuation as seen in the microwave sky 
from a texture collapsing at redshift z ~ 100 a distance r^to — t c away from the observer. 
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